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We construct a map of deflections of ultra-high energy cosmic rays by extragalactic magnetic fields using 
a magneto-hydrodynamical simulation of cosmic structure formation that realistically reproduces the posi- 
tions of known galaxy clusters in the Local Universe. Large deflection angles occur in cluster regions, which 
however cover only an insignificant fraction of the sky. More typical deflections of order ^ 1° are caused by 
crossings of filaments. For protons with energies E > 4 x 10 19 eV, deflections do not exceed a few degrees 
over most of the sky up to a propagation distance of 500 Mpc. Given that the field strength of our simulated 
intergalactic magnetic field forms a plausible upper limit, we conclude that charged particle astronomy is in 
principle possible. 
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ON Introduction. Considerable effort is presently under- 
.taken around the world to build ^ |2] experiments de- 
'voted to determining the composition, the energy spec- 
.trum and the arrival directions of Ultra High Energy 
jCosmic Rays (UHECR). This challenge is in part mo- 
tivated by the Greisen-Zatscpin-Kuzmin (GZK) puz- 
zle [3] which became particularly acute with Fly-Eye 
■and AGASA data [2], and by the realization that the 
UHECR flux at E > 10 19 eV is probably dominated by 
the emission of sources which are quite different from 
^ .conventional galactic sources. The directional informa- 
k™3 tion may allow the identification of UHECR sources, 
5_j 'provided primary particles are not deflected too much 
'by galactic and intergalactic magnetic fields (IGMFs). 

Several arguments suggest that UHECR are electri- 
cally charged nuclei, most probably they are protons. 
The possibility of neutral particles is not ruled out, but 
needs not be discussed here since such rays point back 
to the sources anyway. It is possible that a fraction of 
UHECR is comprised of iron nuclei, see e.g. Ref. [S]. 
However, according to an analysis of inclined events 
recorded by the Haverah Park shower detector |S] , above 
10 19 eV less than 30% of the primary cosmic rays can be 
iron nuclei at the 95% confidence level. In what follows 
we normalize our results to the case of protons. The 
case of other nuclei can be recovered by multiplication 
with their charge, Z. 

Galactic magnetic fields (MF) with B ga \ ~ 1 fj,G are 
not expected to produce significant deflections at ex- 
tremely high energies, E > 10 20 eV in the case of pro- 
tons. Even at lower energies i?~4x 10 19 eV, strategies 



have been proposed which allow source identification 
without detailed knowledge of the galactic MF [7|. 

The very attractive perspective to do astronomy with 
proton primaries might however be spoiled by the pres- 
ence of strong IGMFs. So far, evidences of the presence 
of IGMFs have been found only within, or very close 
to, rich clusters of galaxies. The most relevant observa- 
tions are those based on Faraday rotation measurements 
(RM) of the polarized radio emission of sources located 
within or behind clusters, and on the synchrotron emis- 
sion of relativistic electrons in the intracluster MF. The 
results of both methods imply the presence of MFs with 
strength at the nG level extending up to 1 Mpc from 
cluster centers. The coherence length of the field is 
inferred to lie in the range 10 — 100 kpc (see recent 
review |S] and references therein). Such fields do cer- 
tainly induce large deflections of UHECR protons that 
cross clusters of galaxies. However, galaxy clusters fill 
only a tiny fraction of the volume of the universe, so 
that we may expect them to produce large deflections 
at best over a small portion of the sky E3]- Out- 
side clusters, only upper limits on the IGMF strength 
are available. They are at the level of 1 — 10 nG for 
fields extending over cosmological distances with coher- 
ence lengths in the range 50 to 1 Mpc, respectively 11 . 
These bounds do not hold for MFs in clustered regions, 
like filaments connecting galaxy clusters where the field 
might be as large as 10 -7 G. In principle, either a weak 
all pervading smooth field, or stronger fields localized in 
a complex web of filaments, may produce sizable deflec- 
tions of UHECR over a large portion of the sky. It is 
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hence evident that a better knowledge of the large-scale 
magnetic structure of the universe is called for. 

In this letter, we approach this problem by means 
of numerical simulations of cosmic structure formation, 
where we combine the collisionless dynamics of the 
dark matter component with the magnetohydrodynam- 
ics (MHD) of the magnetized gas. Our basic hypothesis 
is that the MFs observed in rich clusters of galaxies are 
the outcome of a MHD amplification process powered 
by the hierarchical formation of clusters. This assump- 
tion is supported by the results of previous simulations 
which, under the same hypothesis that we adopt here, 
succeeded to reproduce the general features of RM in 
several observed clusters ^21- The tiny seed field re- 
quired to initiate the amplification process may be ei- 
ther of primordial origin |13j or the result of a battery 
associated with the initial stages of structure formation 
d 

Simulations of the magnetic structure of the universe 
and of the UHECR propagation within it have been 
previously attempted by several authors ^] 1151 ITB| . A 
novel achievement of our work is that we have performed 
constrained simulations, which reproduce the observed 
large scale structures in the nearby universe, leaving 
essentially no ambiguity for the choice of observer posi- 
tion. This is quite relevant in the present context, since 
it has been shown [^SJ] that the angular distribution, 
as well as the energy spectrum, of UHECRs reaching 
an observer located in a weakly magnetized region may 
differ considerably from that seen by a strongly mag- 
netized one. Furthermore, by tracing UHECR trajec- 
tories in the simulated magnetic structures we are able 
to construct maps of expected UHECR deflections as a 
function of distance that, for the first time, account for 
the actual large-scale structure as seen from the Galaxy. 

MHD simulations of the Local Universe. We use 

initial conditions that were constructed from the IRAS 
1.2-Jy galaxy survey by first smoothing the observed 
galaxy density field on a scale of 7 Mpc, evolving it 
linearly back in time, and then using it as a Gaussian 
constraint for an otherwise random realization of the 
ACDM cosmology. In [Fi^ it was shown that these con- 
strained initial conditions, when evolved to the present 
time, reproduce the observed density and velocity field 
of the Local Universe. In addition, they allow a direct 
identification of prominent clusters (Virgo, Coma, Cen- 
taurus, Hydra, Perseus, A3627, and Pavo) with coun- 
terparts formed in the simulation, which are found at 
the right places, and with approximately the correct ob- 
served masses. We extended the initial conditions of ^Hj 
by adding gas, together with an initial MF. The volume 



filled by high resolution particles within our simulation 
is a sphere of radius ~ 115 Mpc centered on the Milky 
Way. This region comfortably includes the entire Local 
Super Cluster (LSC) and is modelled with a maximum 
spatial resolution of 10 kpc. The simulation uses 51 mil- 
lion gas particles of mass 6.9 x 10 s M , 51 million high- 
resolution dark matter particles, and an additional 7 
million boundary particles in the distant low-resolution 
region. 

We evolved the initial conditions with the newest 
version of the GADGET-code 17:, adding the Mag- 
netic Smoothed Particle Hydrodynamics (MSPH) tech- 
nique 52 to follow MF evolution. Previous work 
[T2"] showed that magnetic seed fields in the range of 
(1 — 5) x 10~ 9 G at redshift z* ~ 20 will be ampli- 
fied due to the structure formation process and repro- 
duce RM in clusters of galaxies. This corresponds to 
B(z*)(l + z„)- 2 ~ 0.2 - 1 x lO" 11 G at the present 
time in the unclustered intergalactic medium (IGM) . It 
was also demonstrated that the MF amplification pro- 
cess completely erases any memory of the initial field 
configuration in high density regions like galaxy clus- 
ters. Therefore, we can safely set the coherence length 
lc(zi n ) of the initial seed field to be infinite in our simula- 
tion. Although this assumption is probably unrealistic, 
it does not lead to underestimation of the UHECR de- 
flections. Concerning the initial strength of the MF, 
we used the highest value which still allowed previ- 
ous MSPH simulations to successfully reproduce RM in 
clusters, i.e. the results presented here give safe upper 
bounds on UHECR deflections. 

Clusters are generally connected by magnetized fila- 
mentary structures of gas and dark matter, where high- 
density filaments often harbor small clusters or groups. 
We find that shock fronts and shear flows are ubiqui- 
tous in these filaments, giving rise to substantial MHD 
amplification in these structures, boosting the MF in- 
tensity far above the expectation of adiabatic compres- 
sion alone, as pointed out in previous work We also 
identified low density filaments where the MF is roughly 
aligned along their axis, with a strength of ~ 10~ 4 /xG. 
This is consistent with a purely adiabatic amplification 
of the seed MF due to the compression of field lines. 
We find no significant MF in the neighborhood of the 
Milky Way's position, which is contrained to lie within 
a sphere of 7 Mpc around the origin. We find a group 
of four halos aligned within the super galactic plane in 
this region, corresponding to the Local Group. Due to 
the lack of small-scale constraints it is not certain which 
of the four galaxies should be best associated with the 
Milky Way, but this does not affect our results. Because 
it is a small galaxy group, the MFs associated with the 
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FIG. 1. Full sky map (area preserving projection) of deflection angles for UHECRs with energy 4 x 10 19 eV using a 
linear color scale. All structure within a radius of 107 Mpc around the position of the Galaxy was used. The coordinate 
system is galactic, with the galactic anti-center in the middle of the map. Positions of identified clusters are marked 
using the locations of the corresponding halos in the simulation. Note that deflections internal to the Milky Way have 
not been included. 



group are sufficiently weak to not lead to significant de- 
flections, despite covering a large fraction of the sky. 

Deflections of charged UHECR. Having obtained 
a 3D model of MFs in the Local Universe, we can con- 
struct an associated map of deflections of charged par- 
ticles under the action of the Lorentz force. We here 
consider only protons with energy E = 4 x 10 19 eV. 
This is the threshold value for the process of photo-pion 
production in collisions with Cosmic Microwave Back- 
ground (CMB) photons (p + "fcMB — > p(n) + 7r°( + )). 
The energy loss length is large, Ie ~ 1000 Mpc (for a 
recent review see e.g. an( i initially higher proton 

energies quickly degrade into this range. Neglecting en- 
ergy losses and taking r? = 4x 10 19 eV to be the energy 
at detection, we obtain upper bounds for the deflections 
of protons with higher energy since the deflection angle 
decreases linearly increasing the energy. 

We do not follow particle trajectories directly; in- 
stead we compute accumulated deflections along recti- 
linear paths. This is a reasonable simplification since 
we are not interested in actual source positions, but 
rather in finding directions with small deflections. In 
Fig. ^ we show a deflections map obtained by tracing 



an isotropic distribution of protons from a maximal dis- 
tance of (imax = 107 Mpc to the observer. Recall that 
Fig. ^represents a map of deflections, not a distribution 
of arrival directions. The former is independent of the 
assumed distribution of UHECR sources. 

The pattern of clusters and filaments is clearly vis- 
ible in Fig. ^ Large deflections are produced only 
when protons cross the central regions of galaxy clus- 
ters, and most of these strongest deflections are found 
along a strip which can be approximately identified with 
the Great Attractor. The observed positions of Virgo, 
Coma, Hydra and Centaurus lie in this region. Their lo- 
cations quite precisely coincide with regions where the 
deflections exceed 4°. Perseus and other minor clusters 
produce large deflections in other well delineated regions 
of the sky. Outside clusters, which occupy only a small 
fraction of the sky, deflections of 1° — 2° occur along an 
intricate network of filaments, covering a larger area. 
The regions with 6 <C 1° correspond to voids where the 
MF strength is even smaller than 10~ n G. 

In order to investigate the relative importance of de- 
flectors at different distance, we also produced deflection 
maps that only included deflectors up to some maximum 
distance. We observe no significant deflections produced 
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at distances smaller than 7 Mpc. Massive clusters at 
large distances (~ 100 Mpc) produce large deflections 
but cover only a negligible fraction of the sky, so that the 
bulk of the deflections is produced by passages through 
filaments. 

In Fig. EJ we plot the fraction of the sky, A(Sth), over 
which deflections larger than <5 t h are found, for different 
propagation distances. We see that deflections larger 
than 1° are to be expected over less than 20% of the 
sky up to the distance d = 107 Mpc. For large dis- 
tances d, we find that ^4(<5th, d) approaches a self-similar 
behavior, viz. ^4(£ t h,d) = A n (5 t h x (do/d) a ). Numeri- 
cally, we observe a = 0.8 for 70 < d/Mpc < 110. Self- 
similarity is consistent with the assumption that the 
density of deflectors (filaments) reaches a constant value 
at large distances. Since MFs are uncorrelated in dif- 
ferent filaments, multiple filament crossings should pro- 
duce a "random walk" in the deflection angle, resulting 
in a = 0.5. The value of a = 0.8 we observe may hence 
indicate that the regime of multiple filament crossings 
is not yet reached over the distances probed by our sim- 
ulation. We include an extrapolation of A(5th, d) up to 
a distance of 500 Mpc in Fig. [21 shown for two values of 
a, the observed one of a = 0.8, and the expected one for 
large propagation distances, a — 0.5. We expect that 
these two curves bracket the range of true deflections at 
E~4x 10 19 eV. 

We comment finally on the potential effect of the un- 
clustered component of the IGMF, i.e. the field in voids 
and low density regions outside of clusters and filaments, 
coherence length of which, / c , is unconstrained by our 
simulations. If l c <C d the proton trajectory makes a 
random walk through the magnetic domains, and the 
overall deflection is given by 
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Hence, observable deflections are not produced by the 
unclustered component of the IGM if l c is smaller than 
a few tens of Mpc. Note that such small coherence 
lengths are expected from most of the proposed gen- 
eration mechanisms of seed IGMFs ^31 ■ The few mech- 
anisms predicting larger l c generally give rise to MF 
which are too weak to produce observable deflections of 
UHECR. Furthermore, an unclustered IGMF does not 
exist at all if the seed field is generated by a battery 
powered by structure formation [14] . 

Conclusions. We presented the first map of UHECR 
deflections in the Local Universe that is based on a sim- 
ulation that realistically reproduces the known large- 
scale structure around the Galaxy, while simultaneously 
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FIG. 2. Cumulative fraction of the sky with deflection 
angle larger than <5 t h, for several values of propagation 
distance (solid lines). We also include an extrapola- 
tion to 500 Mpc, assuming self similarity with a — 0.5 
(dashed line) or a — 0.8 (dotted line). The assumed 
UHECR energy for all lines is 4.0 x 10 19 eV. 



following the MHD amplification of MFs during cos- 
mic structure formation. The positions and masses of 
the most prominent clusters are reproduced well in our 
simulation. This is an important advantage of our tech- 
nique. Since local structures subtend large angles on the 
sky, it is important to be able to reliably identify "bad" 
regions of expected large deflections, a task that can be 
accomplished using our map, thereby providing impor- 
tant guidance for UHECR source identification. Pro- 
vided our basic hypothesis about the origin of IGMF 
is correct, our results should be understood as upper 
bounds for the expected deflection angles, because we 
have used the largest seed field still compatible with 
the RM in clusters, and secondly, we neglected UHECR 
energy losses on the path to the detector. The actual 
observation of stronger deflections would imply that the 
evolution of the IGMF is not merely passive, possibly 
indicating a pollution of the IGM by physical process 
such as galactic winds. 

We have also extrapolated the distribution of deflec- 
tion angles to very large source distances in a statistical 
manner. Out to 500 Mpc and at E > 4 x 10 19 eV, 
typical deflections are smaller than the angular reso- 
lution of current ground array UHECR detectors over 
more than half of the sky (but may exceed angular 
resolution of stereo fluorescent detectors). This result 
is consistent with an observed small-scale clustering of 
UHECR arrival directions |3U] ^ and with evidence for a 
BL Lacs - UHECR correlation |23] in the energy range 
E ~ 4 x 10 19 eV being due to protons j2S|. On the 

There is no clustering in the current HiRes data I21II22I . which 
became available after submission of our paper. However, with the 
current statistics there is no contradiction yet 1221 1231 . 
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other hand, our results do not support models which 
invoke strong MFs in the local universe to solve the 
GZK anomaly as well as models which explain small- 
scale clustering by magnetic lensing. 

We conclude that charged particle astronomy should 
in principle be possible regardless of the way the GZK 
problem will be resolved. 
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